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ABSTRACT 

Young, massive stars dominate the rest-frame ultraviolet spectra of star-forming galax- 
ies. At high redshifts (z > 2), these rest-UV features are shifted into the observed- 
frame optical and a combination of gravitational lensing, deep spectroscopy and spec- 
tral stacking analysis allows the stellar population characteristics of these sources to 
be investigated. We use our stellar population synthesis code BPASS to fit two strong 
rest-UV spectral features in published Lyman-break galaxy spectra, taking into ac- 
count the effects of binary evolution on the stellar spectrum. In particular, we con- 
sider the effects of quasi-homogeneous evolution (arising from the rotational mixing 
of rapidly-rotating stars), metallicity and the relative abundance of carbon and oxy- 
gen on the observed strengths of He II A 1640A and C IV A 1548,155lA spectral lines. 
We find that Lyman-break galaxy spectra at z ~ 2 — 3 are best fit with moderately 
sub-solar metallicities, and with a depleted carbon-to-oxygen ratio. We also find that 
the spectra of the lowest metallicity sources are best fit with model spectra in which 
the Hell emission line is boosted by the inclusion of the effect of massive stars be- 
ing spun-up during binary mass-transfer so these rapidly-rotating stars experiencing 
quasi-homogeneous evolution. 

Key words: galaxies: starburst - galaxies: stellar content - binaries: general - stars: 
evolution - stars: Wolf-Rayet 



1 INTRODUCTION 

Massive stars dominate the light and kinematics of young 
and star-forming galaxies. They spend most of their life- 
times with effective temperatures greater than 10 4 K, and 
hence most of their light is emitted in the ultraviolet (UV) 
region of the electromagnetic spectrum. It is difficult to ob- 
serve emission in the far-UV (~1000-2000 A) from sources 
in the local Universe, due to absorption in the Earth's atmo- 
sphere and the limitations of space-based instruments sen- 
sitive to this wavelength range. At redshifts above two this 
emission moves into the range of optical spectroscopy. The 
study of rest-frame UV spectral lines therefore become rela- 
tively straightforward and provides a direct insight into the 
massive stellar population of star-forming and ultraviolet- 
luminous galaxies. 

However a problem arises from the extreme distance 
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and consequent apparent faintness of galaxies at such red- 
shifts. The low signal-to-noise obtained in observations of 
typical individual galaxies makes their line strength and stel- 
lar population properties difficult to determine with any de- 
gree of precision. One approach to circumventing this prob- 
lem and characterising a galaxy population is to combine a 
number of spectra to obtain a single, mean composite spec- 
trum. Studies of su ch stacked spe c tra, n otably the compos- 
ite constructed by IShaplev et al.l (2003) of z ~ 3 rest-UV 
selected 'Lyman-break galaxies' (LBGs), allow the determi- 
nation of absorption and emission line intensities that are 
not detectable in individual sources, while being limited to 
yielding information on a 'typical' source rather than any 
individual galaxy. 



However, the analysis of IShaplev et al.l (|2003T ) revealed 
problems with our understanding of the stellar populations 
in such systems. The s tandard stellar popula tion models pro- 
duced by Starburst99 jLeitherer et al.|[l999l i were not able to 
simultaneously reproduce the observed CIV A1548,155lA 
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stellar P-Cygni absorption feature and Hell A 1640A stel- 
lar emission line. The shape and strength of these broad 
lines are determined primarily by the stellar content of a 
galaxy, in combination with a relatively narrow nebular 
component. Massive main-sequence OB stars dominate the 
C IV line (with a small contribution from highly stripped 
WC Wolf-Rayet stars), while the Hell line emission primar- 
ily arises from evolved, hydrogen poor, Wolf-Rayet (WR) 
stars (with a sm all contribution from rare but luminou s Of 
stars). Recently iBrinchmann. Pettini fc Chariot] l|2008h re- 
vised predictions for the Hell line strengths expected from 
a stellar population by including updated Wolf-Rayet star 
line equivalent widths (EWs) in the Starburst99 population 
synthesis code. From this the authors were able to predict a 
Hell EW tha t was in agreement with that observed in the 
IShaplev et all (|2003l ) spectrum, suggesting that careful con- 
sideration of these massive stars and their evolutionary path- 
ways is vital in the interpretation of these distant sources. 

While the use of composite spectra is necessary to 
study the majority of distant sources, there are now a small 
but growing number of individual Lyman-break galaxies for 
which high signal-to-noise rest-UV spectra have been ob- 
tained. The majority of these observations are made possible 
by high magnification due to strong lensing by intervening 
material along the line of sight. They have highlighted a 
diversity in the population that is obscured by use of a com- 
posite, with some galaxies having extremely strong Hell 
emission and others vanishingly weak, undetectable lines. 
This variety leads to a dilemma in its interpretation. While 
strong C IV lines indicate the presence of massive O stars, 
a weak Hell line would suggest a surprising absence of WR 
stars in the same galaxies. This dichotomy suggests that a 
reexamination of the spectral synthesis predictions of these 
C IV and He II features is necessary. 

Much work has been done over the past few years im- 
proving stellar population and spectral synt hesis and quan- 
tifyin g the uncertainties of t hese m o dels (e.g. Leitherer et al.l 
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approach has been to consider the uncertainties inherent in 
the use of only single-star evolution models as the build- 
ing blocks of synthesis codes. We have expanded our in- 
put stellar evolution models to take account of binary stars 
and the many new evoluti onary paths ways this entails. In 
lEldridge fc Stanwavl (|2009l ) we showed that when recent sin- 
gle star stellar models are used to predict the Hell line 
strength, the measured line equivalent width (EW) is very 
small, and that binary evolution models are a better match 
for the observed Hell strength (and for other spectral WR 
population indicators) in star-forming galaxies. 

Here we build on that work, and consider two uncer- 
tainties in the s tellar evolution mode ls in an effort to in- 
terpret both the IShaplev et al.l (|2003l ) composite spectrum 
and those of individual lensed LBGs at z = 2 — 3. First, we 
investigate the effect of variation in the surface car bon abun- 
dance of OB stars on the derived stellar spectra. lErb et all 
(|201(t) show that the amount of carbon in high redshift 
galaxies, as measured by nebular emission lines, decreases 
more rapidly than their [0/H]-derived metallicity at low 
abundances. Therefore we investigate the effect of decreas- 



ing the relative abundance of carbon on the C IV absorption 
line, in order to examine and model this trend. 

Second, we further investigate the importance of du- 
plicity in a synthesised stellar population when determining 
the He II line equivalent width, assuming both instantaneous 
burst and constant star-formation histories. We now include 
a new evolutionary process that may only be possible at 
low metallicities. Quasi-homogeneous evolution (QHE) oc- 
curs when a massive star rotates sufficiently rapidly that 
the s tar becomes fully mixed during its main-sequence life- 
time (|Maedeijll987l ; lYoon fc Langerl2005l ; iMevnet fc Maederi 
120071 ). While it is difficult for a single star to be born 
rapidly rotating, efficient mass-transfer in massive binaries 
can spin up secondary stars to high rotation rates with ease 
l|Cantiello et al.ll2007l ). At Solar metallicities such stars are 
likely to spin down quickly due to strong stellar winds. How- 
ever with the weakening of stellar winds at lower metallici- 
ties, efficient rotational mixing can occur. In this paper we 
model the effect of this unusual form of evolution on stars 
that accrete material during a mass-transfer event, and in- 
vestigate the dramatic effect such stars would have on the 
integrated spectrum of a stellar population. 

This paper is organised as follows. In Section [2] we dis- 
cuss the modifications made to our Binary Population and 
Spectral Synthesis (BPASS) code in order to investigate the 
phenomena outlined above. In Section [3] we describe the ef- 
fects of these phenomena on the strength of the diagnostic 
He II and C IV spectral features. In Section [4] we go on to 
compare our predicted spectra to a composite spectrum con- 
struct ed from ~900 Lyman-break galaxies bv lShaplev et al] 
(2003) and to various individual examples of z ~ 2 — 3 
Lyman-break galaxies. Finally, in Section [5] we briefly dis- 
cuss our results and outline our conclusions. 



2 SYNTHETIC SPECTRA OF STELLAR 
POPULATIONS 

2.1 Binary Population and Spectral Synthesis 
(BPASS) 

The synthetic spectra used in this paper are cre- 
ated using the Binary Population and Spectral Syn- 
thesis (BPASS) codcJ . It i s described in detail in 
lEldridge. Izzard fc Tout! <l2008l) lEldridge fc Stanwavl (|2009l ) 
and lEldridge. Langer. fc Tout! <|201ll ). 

We use stel la r mod el s from the Cambridge S TARS 
code (|Eggletonl Il97ll ; lEldridge. Izzard fc Tout) 120081 . 
and references therein) , speci fically those calculated in 
lEldridge. Izzard fc Tout) (120081 ). Their key feature is an 
extensive set of detailed binary star models (in addition 
to the detailed single star models) which are key to pro- 
ducing a realistic synthetic stellar population. We consider 
stellar models at five different metallicities: Z = 0.001, 
0.004, 0.008, 0.020 and 0.040 (where a metallicity of 
Z = 0.020 is conventionally considered Solar), with hydro- 
gen mass fraction, X = 0.75 — 2.5Z, helium mass fraction, 
Y — 0.25 + 1.5Z and a default metal distribution given by 
scaled-Solar abundances. 

Given that stellar evolution is non-linear and binary 
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evolution is even less predictable, we do not interpolate be- 
tween models with different masses and initial binary param- 
eters, but rather weight each stellar model by an initial mass 
function (IMF) and distribution of binary properties. Binary 
population fractions, p athways and mass loss ra t es are cal- 
culated as described i nlEldridge. Izzard fc Tout! HH) and 
lEldridge fc Stanwavl |2009l ). The details of our binary in- 
teraction algorithm are relativ ely simple compared to the 
scheme outlined in, for example. iHurlev. Tout fc Pols! (|2002l ) 
and interactions depend primarily on the initial separation 
and stellar mass ratio. Our aim was to investigate the effect 
of enhanced mass loss due to binary interactions on stel- 
lar lifetimes and populations; therefore we concentrated on 
these aspects rather than incorporating additional physical 
processes, each of which would add more free parameters to 
our models and potentially associated uncertainties on those 
parameters or the mechanisms concerned. 



2.2 Quasi- homogeneous evolution 

We include one new evolutionary path, not in the standard 
picture of binary evolution for our secondary stars, first in- 
trodu ced into our population synthesis and described in de- 
tail in lEldridge. Langer. fc Toutl (|201lh . This path only oc- 
curs if the more massive primary star in a binary overfills 
is Roche Lobe and transfers mass to the secondary star. If 
this happens and the secondary has a metallicity mass frac- 
tion less than or equal to 0.004, a mass after accretion of 
more than 10 Mq and has accreted more than 5 percent of 
its initial mass, we assu me it evolves fully mixed during its 
main - sequence lifetime ([Petrovic. Langer fc van der HuchtJ 
2005; ICantiello etail l2007h . This is referred to as quasi- 
homogeneous evolution and is the result of rapid rotation 
due to the a ccretion of material fr om the p rimary star as 



described by Yoon fc Langerl (l2005h (see also Maederl 



Mevnet fc Maederl 120071 ; lYoon. Langer. fc Normanl 
Cantiello et al.ll2007f ). ~ 
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In this case we use simple models in which we assume 
that stars meeting the criteria for QHE are fully mixed dur- 
ing their hydrogen burning evolution. They burn all their 
hydrogen to helium and increase in surface temperature as 
the burning progressing, evolving away from the zero-age 
main-sequence in the wrong direction. Once hydrogen burn- 
ing ends this mixing also comes to an end. This evolution 
prescription is applied whether or not the binary remains 
bound after the primary supernova. The stars never grow to 
large radii, because the homogeneous evolution means the 
increasing molecular weight of the stellar material makes 
the stars shrink throughout their main-sequence evolution. 
Therefore QHE stars in a binary never fill their Roche lobe 
and interact with the other binary component. Other than 
this inclusion of mixing, the stars are treated identically to 
the non-QHE models with the same mass-loss rates applied. 



2.3 Accounting for relative carbon abundance 

In our standard models, and those of other population syn- 
thesis codes, Solar-scaled relative abundances are normally 
employed for a range metals. Variation in the relative abun- 
dance of different species is also possible - for example, an 
enhanced abundance of elements generated by the alpha pro- 



cess is characteristic of bursty, massive star formation. Re- 
cent studies modelling the interstellar medium in some high 
redshift galaxies have suggested that the important ratio of 
carbon to oxygen abundance decreases along with the total 
metallicity of the galaxies fsee lErb et al.|[201ol . for example). 
These abundance variations must be considered in order to 
develop an accurate model for stellar population synthesis. 
However, adapting the synthetic stellar populations is far 
from straightforward. For a full analysis, the stellar models 
and atmosphere models should all be recreated with freely- 
varying relative abundances. This would be a challenging 
and time consuming task, as well as providing many addi- 
tional free parameters which are likely to be unconstrained 
by the data. In this paper we outline a simple prescription 
for estimating the effect of varying carbon/oxygen ratios on 
the synthetic spectrum predicted by a stellar population in 
the UV. 

In addition to variations in the initial carbon/oxygen ra- 
tios, the relative abundance of different elements in a stellar 
interior changes over time due to nuclear reactions. While 
these reactions occur at the core of the star, various pro- 
cesses can make nuclear-processed material observable at 
the surface of a star. For stars above a few Solar masses 
main-sequence hydrogen burning is catalysed by CNO ele- 
ments. During this catalysis the rate determining step is a 
decay of 14 N and so most of the initial CNO abundance ends 
up in this element. This abundance change is concentrated 
at the core of the star, but mass-loss of the outer layers 
and additional mixing (beyond that predicted by mixing- 
length theory) can draw the the processed core material to 
the surface thus changing the observable abundance. Typ- 
ically this is observed as nitrogen enrichment and used to 
infer the importance of ro tational mixing in stellar interiors 
jHunter et al.ll2007l . 120091 ') . 

To take account of a possible decrease in ini- 
tial abundances we have calcul ated WM-Basic 
(|Pauldrach. Hoffmann, fc Lennonl 1200 ll ) models of O 
stars and varied the carbon abundance of the models, 
decreasing the mass fraction of carbon by factors of 0.5, 
0.1 and 0.01. We then measure the equi valent width of 
the u ltraviolet C IV line as described by ICrowther et al.l 
(2006), and compare it to that measured at a standard 
Solar-scaled carbon abundance. The results are presented 
in Table [1] It is apparent that there is not a linear trend 
with metallicity, but rather that the greatest effect is seen 
when Z — 0.004 (i.e. significantly sub-Solar, but not at the 
lowest metallicity). Here even a small decrease in the carbon 
abundance decreases the C IV EW by a substantial fraction. 
At higher metallicities much more carbon must be removed 
before any great effect is seen. At the lowest metallicities 
there is already too little carbon in the atmosphere to 
produce a strong C IV line and so any further reduction 
has a minimal effect. 

To include the effect of varying composition in our spec- 
tral synthesis we use the results in Table [1] to reduce the 
C IV EW in each model atmosphere according to its carbon 
abundance. We found the accounting for the decrease of the 
carbon abundance predicted by stellar models gave rise to 
only a tiny change in the C IV EW. Therefore we applied the 
maximum decre ase in the carbon to oxygen ratio found by 
lErb et al.l l|2010l ) at z ~ 2 so as to investigate the variation 
of the initial carbon/oxygen abundance ratio. We present 
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Table 1. The relative decrease in the CIV line flux at 1550 A 
when the carbon abundance is decreased by a factors of 2, 10 
and 100, The values were estimated from WM-BASIC mod- 
els ijPauldrach. Hoffmann, fc Lennonll200lh with varying carbon 
abundances and constant overall metallicity. X c is the initial car- 
bon mass fraction at each metallicity. 
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Figure 1. The predicted equivalent width of CIV at different 
metallicities measured by [O/H]. Lines represent BPASS models 
assuming constant star formation with the vertical extent du e 
to age. The points are observations from lCrowther et al. I d2006f) . 
Red lines are for single stars, and yellow for binary models, the 
green lines indicate binary models including QHE and the blue 
lines are for binary models with the C IV abundance decreased 
by a factor of four. 



two sets of models, one with normal Solar-scaled carbon 
abundances and a second with the carbon abundance de- 
creased by a factor of four. We only perform this decrease for 
main-sequ ence OB stars when we use the O B model atmo- 
spheres of ISmith. Norris""fc Crowtherl (|2002D . The result is 
general reduction in the strength of the C IV absorption line 
as measured by the equivalent width. As Figure [T] demon- 
strates this reduction can go to explain some of the scatter 
in the observed EWs of C IV in nearby galaxies, although 
the majority of the scatter is due to variations in their star 
formation history. 



2.4 Producing a total synthetic population 
spectrum 

The procedure outlined above yields a synthetic spectrum 
appropriate to each time-step of a stellar evolution model. 
We can then combine the spectra for each star together to 
produce the integrated spectrum for a synthetic stellar pop- 
ul ation. To do thi s we use the initial mass function described 
bv lKroupal |2002i ). This uses an IMF power-law slope of -1.3 




Figure 2. How the strength of the Hell and CIV line vary with 
age and metallicity for various BPASS stellar population models. 
The solid lines show binary models, the dash-dotted lines are 
for binary models including QHE and the dashed lines for single 
stars. Black indicates models with Z = 0.001, the red lines have 
Z = 0.004, the yellow line are for Z = 0.008, the green line for 
Z = 0.020 and the blue line for Z = 0.040. On the binary model 
tracks, plus signs indicate the values at 1 Myr, with crosses along 
the line indicating ages of 10 Myr and 100 Myr , then the asterisk 
indicating the age of 1 Gyr. The light blue boxes indicate the EWs 
measured from the observed spectra described in Section 14.11 



between 0.1 and 0.5 M©, and a slope of -2.35 from 0.5 to 
120 M . 

Finally in our spectral synthesis we include the con- 
tribution from nebular emission. In star-forming galax- 
ies, interstellar gas is ionised by the stellar contin- 
uum emitted blueward of 912A, and upon recombina- 
tion it emits a nebular continuum. Neglecting this emis- 
sion would lead to an incorrect estimate of the equivalent 
widths of emission lines and incorrect broad-band c olours 
jZackrisson, Bergvall fc LeitetJliool ; iMolla et al.ll2009l). We 
use t he radiative transfer program Cloudy ( Ferland et al.l 
1998) to produce a detailed model of the output nebular 
emission spectrum excited by our stellar spectra. The model 
output is sensitive to the chosen geometry, inner radius and 
composition of the gas used in the code. The details of our 
illustrative nebular emissio n model are identical to those in 
lEldridee fc Stanwavl l|2009l ), and we output the final contin- 
uum and line strengths for use in our synthetic spectra. 



3 PREDICTED EQUIVALENT WIDTHS 

We consider in this work the CIV and Hell spectral lines, 
both of which exhibit a broad component whose strength is 
determined by the stellar spectrum. We have already shown 
in Figure [1] that model predictions for the C IV line agree 
with those deduced from observations. Conventional wisdom 
has been that the magnitude of He II line EW decreases 
in tandem with the equivalent width of C IV as shown by 
iBrinchmann. Pettini fc Chari ot (2008). This is because it is 
assumed that at lower metallicity there are fewer WR stars 
with shorter lifetimes to contribute to the spectral features. 
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Figure 3. The predicted equivalent widths for the CIV and Hell stellar wind lines (upper and lower panels respectively). The panels in 
the left-hand column assume an instantaneous burst of star formation. The panels in the right column assume a constant star formation 
rate. The solid lines are for the standard binary models, the dashed line are those from single stars models and the dash-dotted lines for 
binary models including QHE. The black lines are for Z = 0.001, the red lines for Z = 0.004, the yellow line for Z = 0.008, the green 
line for Z = 0.020 and the blue line for Z = 0.040. 



However, as shown in lEldridee fc Stanwavl l|2009h . the inclu- 
sion of binaries make this predicted relation less clear-cut. 

In Figure [2] we show how the C IV EW, which does 
decrease with metallicity, relates to the Hell line strength 
for our different model populations. When only single star 
populations are considered, the strength of the Hell line is 
largely independent of C IV line strength, and, surprisingly, 
peaks at around Z — 0.008 rather than at the highest (Solar 
and super-Solar) metallicity. This is because absorption lines 
in B star spectra decrease the apparent strength of Hell 
at higher metallicity. Stellar populations including binaries 
reach higher EWs than the single star populations because 
of their greater number of WR stars (produced due to mass- 
transfer events in the binary systems). The highest EWs are 
only possible, however, when the effect of QHE is included in 
a binary population. The production of long lived hydrogen 
burning stars that are observed as WR stars boosts the line 
strength at the lowest metallicities. 

The same trends are evident in Figure [3] The highest 
Hell EWs (i.e. strongest line emission) arise at later ages 
in QHE models than for the single star populations. This 
is because the additional WR stars produced from binary 



evolution and QHE are typically from lower mass stars that 
would not normally become WR stars and take of order 10 7 
years to evolve to this state. Also in Figure [3] we show how 
both C IV and He II lines vary assuming an instantaneous 
burst of star formation. This model may be more appropri- 
ate at the highest redshifts (z > 5) where the stellar popu- 
lations appear to be short-lived, dramatic eve nts with little 
fuel available for ongoing star formation (see iDavies et al.l 
l20ld ). 



In this instantaneous star formation case, larger EWs 
are possible for Hell than in the constant star formation 
model. This is only because there are no younger stars con- 
tributing to the continuum and reducing the relative sig- 
nificance of WR line flux in the equivalent width (i.e. line- 
to-continuum flux ratio). However the Hell emission lines 
typically observed in Lyman-break galaxies with good spec- 
troscopic data at z ~ 2 — 3 do not require these extreme 
EWs, but remain in agreement with those predicted by our 
constant star formation models. 
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4 APPLICATIONS TO THE REST-FRAME UV 
SPECTRA OF LYMAN-BREAK GALAXIES 

4.1 Observations 

We have collated the spectra of five Lyman Break galaxies 
at 2 = 2 — 3 for which high quality rest-UV spectroscopy 
has been published, four are gravitationally-lensed galaxies 
and the other spectrum is the result of a deep, 12 hour inte- 
gration. The strong magnification of distant sources by the 
gravitational effects of intervening matter allows the study 
of individual galaxies, rather than the necessary blurring 
of any variation in a population through use of a compos- 
ite. Before comparing these spectra to models it is useful to 
consider the spectra in co mparison with each o ther and the 
earlier LBG composite of IShaplev et al.l l|2003h . 

The data, in order of publication, comprise rest- 
frame ultraviolet sp ectra for MS 1512-cB58 (hereafter cB58, 
IPettini et al. 20021) at z = 2.72, the Cosmic Horseshoe 



( Belokurov et al.ll2007i; iQuider et al1l2009l) at z = 2.38, the 



8 o'c lock arc ( Allam et al.ll2007l ; iDessauges-Zavadskv et al.l 
bold ) at z = 2.73, the Cosmic Eve (|Quider et al.ll2010l) at 
z = 3.07 and Q2343-BX418 (hereafter BX418. lErb et all 
l2010l ) at z = 2.30. In Figure]?] we combine the normalised 
rest-frame ultraviolet spectra to form a composite of the 
lensed galaxies and compare it to the z m 3 Lyman-break 
galaxy composite of Sha plev et al.l (|2003T ). As the upper pan- 
els of the figure illustrate, there appears to be a remarkable 
qualitative agreement between our small composite and that 
based on ~ 900 individual galaxies at lower signal to noise. 
This suggests that, taken as a group, these Lyman-break 
galaxy examples are representative of the more numerous 
photometrically-selected sample. The primary difference be- 
tween the two composites is that the C IV absorption feature 
(in both its broad stellar wind and narrow nebular compo- 
nents) is somewhat deeper in the lensed galaxies than, on 
average, in their unlensed counterparts. 

However, even this small composite does not tell the 
whole story. When studying the C IV and He II profiles of 
the individual spectra it is apparent that the lensed sources 
form two loose groups that can be separated by the strength 
of their Hell emission line. One set contains strong Hell 
stellar emission with equivalent wid ths (EWs) of 2A or more 
I Allam et al.ll2007l ; lErb et aljfeoid ). while a second group 
presents only weak or absent Hell emission (|Pettini et al.l 
12002] ; iQuider et al.ll2009l . l2010l ). We compare the mean spec- 
tra of these high and low Hell EW sets in the second set 
of panels in Figure [?] Not only does the He II strength vary 
significantly between these two cases, but the C IV profile 
is also deeper in the weak-Hell case, indicative of a higher 
carbon abundance. 

In the final two sets of panels we show the individual 
spectra of the lensed galaxies. It is clear that the spectra of 
sources in the low He II EW category are remarkably uni- 
form. There is very little difference between any of the ob- 
served spectra. The two sources with high He II EW spectra 
are also similar to one another and differ primarily in the 
strong (relatively narrow) nebular e mission contribut ing to 
the He II line of the BX418 spectrum. lErb et all (|20ld ) fit a 
dual Gaussian model to the He II emission line in this source 
and estimate that 20 per cent of the line flux is contributed 
by nebular emission. Neglecting the nebular contribution, 
the two spectra have identical He II equivalent widths. 
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Figure 4. Comparison of the different observed spectra. In the 
upper two panels we compare the Shapley composite with a com- 
posite of the lensed galaxies, below this we compare composites of 
the high and low Hell galaxies, in the final we panels we compare 
each lensed spectra to the relevant composite. 



This qualitative comparison suggests that there may be 
evidence for two distinct classes of z ~ 3 LBGs, separated 
by the strength of the broad He II stellar emission line. Fur- 
thermore, it would appear likely that both types of galaxies 
are represented within the photometrically-selected galaxy 
sample contributing to the Shapley et al. spectral stack. 



4.2 Modelling observed z ~ 2 — 3 galaxies 

4-2.1 Shapley et al (2003) composite 

Given the new elements in our stellar population synthesis 
code, it is informative to perform a direct comparison be- 
tween the spectral features predicted by the refined BPASS 
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Figure 5. A comparison between spectral featu res in the ob- 
served composite rest-frame UV spectra LBG of IShaplev et all 
ll2003ll in blue and our binary population models including QHE, 
shown with and without a reduction in the carbon abundance 
in red and black respectively. Left-hand panels are for C IV and 
panels on the right for the Hell line. The model spectra are for 
continuous star-formation lasting 10 8 ' 5 yrs. We show four differ- 
ent metallicities, from top to bottom Z = 0.001, 0.004, 0.008 and 
0.020. 



models and those observed in the six Lyman-break galaxy 
spectra (five individual, one composite) described in Section 
14.11 This comprehensive study enables us to gain further in- 
sight into both the accuracy of our models and the implied 
properties of the high-redshift galaxies themselves. 

In Figure [S] we present a comparison between the com- 
posite spectrum of ^900 LBGs with a mean redshift of 
three jShaplev et alj l2003h and our model spectra, taken 
log(age/yrs)=8.5 into an ongoing star formation episode 
(found to be the median age and st ar formation history 
of this sample by IShaplev et all l200lh ■ In the original pa- 



Figure 6. Similar to Figure[5]but showing a pure single-star pop- 
ulation in black and a binary population model excluding QHE 
in red. 



per presenting this rest-UV composite, the authors found 
it difficult to simultaneously reproduce the line strengths 
of the C IV and He I I line without resorting to mode l s with 
extreme IMFs. Later lBrinchmann. Pettini fc Charlotj (|2008h 
was able to predict the required EW for the Hell line but 
only at a lower metallicity than suggested by other measure- 
ments of the sample. 

We see in Figure [S] that the C IV line in the LBG com- 
posite is relative weak compared to our Solar metallicity 
models. Only at sub-Solar metallicities (Z—0. 001-0. 004) can 
the models reproduce the broad component of the observed 
P-Cygni wind profile shape. However at the higher of these 
metallicities, the line profile can only be reproduced if the 
C/O ratio is a quarter of that in the standard Solar abun- 
dances. This is close to the depletion of this a bundance ratio 
by a factor of 0.2 (relative to Solar) found bv lShaplev et all 
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(2003) for the nebular gas in this source. This implies that 
the stars we are observing likely formed from the same gas 
that still surrounds them at the epoch of observation. 

For completeness, we compare the composite spectrum 
to our theoretical single star and non-QHE binary evolution 
models in Figure [6] There is little difference in the C IV line 
relative to the models incorporating QHE shown in Figure 
[5] However both our standard single star and earlier binary 
models have very weak He II lines compared to those spectra 
where we apply a prescription for QHE, and under-predict 
the observed line strength. This indicates that accounting 
for binary stellar evolution alone cannot entirely explain the 
observed Hell profile, but rotation in binaries may play an 
important role in producing the observed spectra. 

Because the spectrum analysed here is a composite of a 
number of different galaxies we can only attempt to deter- 
mine a mean or range of metallicities. A qualitative analysis 
suggests an allowed metallicity range of 0.001 < Z < 0.008 
for the 'typical' Lyman-break galaxy at z ~ 3. This corre- 
sponds to a 7.6 < [O/H] < 8.5 (or 0.1 < Z/Z e < 0.6 where 
[0/H]q = 8.7), and is consistent with the pr eviou s measure- 
ments determined by IShaplev et alj (|2003l ) and iRix et all 
l|2004h . 

4-2.2 Individual galaxies 

We address the spectra of the individual lensed galaxies com- 
piled here starting with the sources with the lowest inferred 
metallicity, based on our model comparisons. In each case, 
we make the comparison at a constant log(age/years)=8.5 
(age=320 Myr), based on the median star- formation derived 
age of the z ~ 3 Lyman-break galaxy population as a whole. 
We note that, as figure [3] indicates, there is little evolution 
in the model strengths of our lines of interest at stellar pop- 
ulation ages above ~100Myr. 

The UV sp ectrum of the z = 2.3 galaxy BX418 from 
lErb et all (|201CT ) (shown in Figure has a very weak C IV 
stellar wind line with a P Cygni profile and strong and broad 
Hell emission. There is clearly a narrow nebular compo- 
nent superposed on the broader Hell stellar emission line. 
This indicates a very hard ionising spectrum present within 
the galaxy. The He II line implies that the metallicity must 
be low, especially in carbon, suggesting a metallicity range 
of 0.001 to 0.004 (i.e. significantly sub-Solar). However the 
more prominent CIV profile implies a slightly higher (but 
still sub-Solar) carbon abundance in this source. As is the 
case for the Shapley et al composite, matching both lines 
simultaneously requires the presence of both a reduced rela- 
tive carbon abundance and QHE in the population synthesis 
model. 

In Figure [5] we consider the 8 o'clock arc, a z — 
2.7 lensed Lyman bre ak galaxy, with a rest-frame UV 
spectrum obtained by iDessauges-Zavadskv et all 42O10l ). 
This galaxy has reported metallicities higher than oth- 
ers in this sample (Solar, or near-Solar) and has also 
been reported as relatively mass ive and dusty (M« ~ 
4 x 10 11 M o , E(B-V)~0 67, s ee iFinkelstein et alj 120091 ; 
IDessauges-Zavadskv et aT1l2010l ). In this context, its rest-UV 
spectrum is somewhat surprising. As in the case of BX418, 
the 8 o'clock arc has a very weak C IV stellar wind line (the 
absorption is dominated by the narrow nebular component), 
but strong and broad Hell emission, without any obvious 



BX418 log(Age/yrs) = 8.5 
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Figure 7. Similar to Figu re \Ej but for the observed UV spectrum 
of BX418 l|Erb et alj|2010h . 

nebular component. The lack of a strong C IV wind line in- 
dicates that the preferred metallicity in the regions dominat- 
ing the rest-UV emission must be low, especially in carbon, 
suggesting a metallicity range of 0.001 to 0.004, with higher 
metallicities within this range requiring a reduced [C/O] ra- 
tio. The strong He II bump suggests metallicities at the high 
end of the allowed range to reproduce this line. At the lower 
metallicities fewer WR stars are produced, and weaker line 
emission in the WR stars also reduces the strength of the line 
observed. This He II line strength is also only reproducible 
with QHE in the stellar population. 

Based solely on the stellar-dominated spectral features, 
the similarity between the rest-UV spectrum of this galaxy 
to that of BX418 (shown in Figure [4j would suggest they 
are of a similar metallicity, with the deeper absorption fea- 
tures hinting at a slightly more metal rich environment in 
the 8 o'clock arc. This is somewhat in contradiction to the 
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Figure 8. Similar to Figurc[5]b ut showing th e observed UV spec- 
trum of the 8 o'clock arc IjAllam et aI.ll2007T> . 



Figure 9. Similar to Figure \Ej but for t he observed UV spectrum 
of the Cosmic Eye l|Ouider et al.l 2010). 



higher deri ved metallicity range found by Finkelstein et alj 
( 2009) and lDessauges-Zavadskv et all (|201o1 T based on nebu- 
lar ionization models and older stellar population models. It 
is possible that any discrepancy may arise in part from the 
different regions probed by these different measurements. If 
so, it would imply that the youngest star formation (and 
hence the massive stars that dominate the He II and C IV 
lines under consideration here) is taking place in regions 
of lower metallicity than the more metal-rich gas elsewhere 
in the galaxy. It is, of course, dangerous to over-interpret 
weak features in data with this signal to noise. However, the 
flattened line profile of the Hell emission may support the 
hypothesis that metal-rich gas is present to some degree in 
this system, if interpreted as the superposition of a narrower 
(perhaps nebular) absorption line on the peak of the broad 

stellar emission fe ature. 

Interestingly, ICresci et al.l (|2010l ) have identified three 



star-forming galaxies at a similar redshift that show evi- 
dence for an inverted metallicity gradient (i.e. central star 
forming regions with lower metallicity than outlying re- 
gions) which may arise from the central accretion of cold 
gas. Like the 8 o'clock arc, all three are relatively massive 
(M* ~ IO^Mq) and they also show dynamics that sug- 
gest rotational support rather than merger-driven star for- 
mation (which may be more dominant in lower-mass UV- 
luminous systems). A cold-accreting system of this kind is 
likely to display the same discrepancy between metallicities 
determined from the youngest stars and from other meth- 
ods. Cold-mode accretion has been hypothesised as a dom- 
inant galaxy form ation mechanism at the highest redshifts 
|Dekel et al.ll2009l ), although th i s hyp othesis remains con- 
troversial (see e.g. ISteidel et al.l[201ol l. While such a case 
has not been proven in the 8 o'clock arc, it may prove an in- 
teresting target for detailed analysis with integral field spec- 
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Figure 10 . Similar to Figur c[5]but for the observed UV spectrum 
of cB58 l|Pettini et al.ll2002t ). 



Figure 11. Similar to Figurc[5] but for the ob served UV spectrum 
of the Cosmic Horseshoe jOuider et al,|[2u09l y 



troscopy, to determine whether there is any evidence for this 
scenario. 

The remaining three galaxies differ to those already dis- 
cussed in that they do not require the simultaneous invoca- 
tion of low metallicity, reduced carbon abundance and QHE 
in their synthetic stellar population to explain the spec- 
tral features under consideration, but rather display some 
or none of these features. 

The z = 3.1 Cosmic Eye spectrum (|Quider et al.ll201(il ) 
is illustrated in Figure [9] and is notably devoid of strong 
He II emission. This can be explained if the galaxy lies above 
the metallicity limit at which QHE no longer occurs. The 
absorption section of the C IV P Cygni profile also implies a 
relatively strong metallicity somewhere between Z = 0.004 
and 0.008, with a reduced carbon abundance required at 
the higher metallicity mass fraction. Higher metallicities are 
ruled out by the relatively shallow absorption. 



The remainin g two examples, cB58 (z = 2.7, 

)2 ) and the Cosmic Horseshoe [z — 2.4, 

iQuider et al.ll2009T ). in Figures llOl and [TT1 respectively, tell 



iPettini et al 



a similar story. The lack of a strong Hell line, yet apparent 
C IV P-Cygni profiles, suggest higher metallicities for these 
galaxies, above the range of metallicities for galaxies with 
stronger Hell. 

These results suggest that the existence of two popula- 
tions of LBG galaxies, defined by the presence or absence of 
broad Hell emission as discussed in Section [4.11 may be ex- 
plained by intrinsic galaxy properties differentiated primar- 
ily by the overall metallicity mass fraction of the youngest 
stellar population (which dominates the UV emission and 
hence the two stellar features under consideration here). If 
the metallicity mass fraction is low then QHE is possible 
in binaries and produces more long lived Of and WR stars, 
which contribute to the broad Hell line. By contrast, when 
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the metallicity is high QHE is no longer possible because 
stellar wind mass-loss breaks the rotation of stars rapidly to 
prevent this unusual evolution, and Hell emission is com- 
paratively weak. Importantly, in both sets of galaxies we 
see evidence for a sub-Solar carbon to oxygen ratio in the 
weaker than expected C IV P-Cygni profiles. 



5 CONCLUSIONS 

In this paper we provide evidence that our new synthetic 
stellar population models can explain the observed stellar 
emission lines in the UV spectra of LBGs, with best fit mod- 
els as summarised in Tabled We see evidence for two main 
types of LBG UV spectra, differentiated by the presence 
or absence of broad Hell line emission, and explained as 
the result of varying metallicity in the youngest stellar pop- 
ulation. This leads to the absence/presence of an extreme 
rotational effect (quasi-homogeneous evolution) in the stel- 
lar binaries incorporated in our synthesis models. We find 
three details that must be included to model the spectra 
accurately: massive binary stars, QHE from the spin-up of 
the secondary stars in binaries in the event of mass-transfer 
events, and a reduced carbon-to-oxygen abundance ratio. 

Our models support the inference that the sub-Solar 
[C/O] ratio derived from n ebula abundances in the z — 2.3 
galaxy studied in detail bv lErb et al.l (|2010h may be wide- 
spread in this population, and that its inclusion in synthe- 
sis models is important when fitting individual spectral fea- 
tures. Depletion by a factor of approximately four is consis- 
tent with the entirety of the sample discussed here. 

In addition, we confirm our earlier work which suggests 
that binary evolution has significant effects on the spectra 
of unresolved stellar populations. The main effect of this is 
increased mass loss in binary stars, and therefore an increase 
in the overall number of WR stars. Furthermore the effects 
of rotation (QHE) within a binary system can also lead to 
dramatic and hither-to unexpected results in the derived 
model stellar spectra. 

Given the near-ubiquity of sub-Solar carbon-to-oxygen 
ratios in this sample, and the increasing import of QHE 
effects at low metallicity, it is interesting to consider the 
possible implications for analysis of Lyman-break galaxies 
at still higher redshift. This analysis has focused on galaxies 
at z < 3. However similar sources, selected for high rest-UV 
flux (and hence recent star formation) have now been iden- 
tified to z — 8 and beyond. The detailed spectroscopic anal- 
ysis of these sources is still more difficult than that at z = 3 
due to a combination of increased sky background at long 
wavelengths, decreased detector sensitivity and increasing 
luminosity distance, and is not considered further here. How- 
ever, forthcoming facilities such as the James Clark Maxwell 
Telescope (JCMT) and the KMOS spectrograph may make 
analysis of the Hell line and other, longer wavelength spec- 
tral signatures relatively straightforward at z = 4 and above. 
At t he low, but non-zero, metallicities infe rred at z — 4 — 8 
fe.g. lAndo et al.ll2004rDouglas et al.ll2010h . the effects of de- 
pleted carbon abundance and QHE may have a significant 
effect on observed spectral features and improved modelling 
in these areas will be essential for accurate interpretation of 
spectra obtained. 

As the case of the 8 o'clock arc makes clear, it is possi- 



ble that the youngest star formation is atypical of the galaxy 
as a whole, which has implications for our understanding of 
the processes of galaxy formation. If cold-mode accretion is 
indeed increasi ngly common at high redshifts as has been 
suggested (e.g. ICresci et af] l20ld ; IPekel et al.1 l2009h then 
stellar emission and absorption features may be an acces- 
sible indicator of this. While the limited sample discussed 
here is, of course, difficult to interpret, the use of KMOS 
and the JCMT to build a sample of high redshift sources 
with high quality rest-UV spectra may allow the identifica- 
tion of galaxies for which the youngest stellar populations 
are anomalously low in metallicity. Such sources may well 
prove interesting targets to study as possible examples of 
the still hotly-debated cold-mode accretion scenario. 
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Table 2. The stellar population model found to best fit each of the Lyman-break galaxies discussed in section ^. 21 indicating the preferred 
range of metallicity and whether either a reduced carbon to oxygen ratio or quasi-homogeneous stellar evolution models are required to 
explain the observed spectrum. We list the details for composite LBG spectrum (which demonstrates the broadest allowed metallicity 
range) and the observed samples with strong and weak Hell line emission as discussed in Section |4, 21 
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